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Abstract 

Three-dimensional Navier-Stokes computational fluid dynamics (CFD) analysis has been 
performed in an effort to determine thermal boundary layer correction factors for circular 
convective heat flux gauges (such as Schmidt-Boelter and plug type)mounted flush in a flat 
plate subjected to a stepwise surface temperature discontinuity. Turbulent flow solutions with 
temperature-dependent properties are obtained for a free stream Reynolds number of 1E6, 
and freestream Mach numbers of 2 and 4. The effect of gauge diameter and the plate surface 
temperature have been investigated. The 3-D CFD results for the heat flux correction factors 
are compared to quasi-2D results deduced from constant property integral solutions and also 
2-D CFD analysis with both constant and variable properties. The role of three- 
dimensionality and of property variations on the heat flux correction factors has been 
demonstrated. 
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Nomenclature 

b (T w \ -Taw)/(T w2 -T aw ) 

c sound speed, *JyRT 

c p specific heat at constant pressure 

h heat transfer coefficient 

k thermal conductivity 

L reference plate length 

M Mach number, /c 

Nu x local Nusselt number, h x x/k 
Pr Prandtl number 

q heat flux, 

R gas constant; also radius of heat flux gauge 

Re Reynolds number 

r recovery factor 

St x local Stanton number, h x l\p„c p u m )= Nu x /(Re x Pr) 

T temperature 

k, v, w velocities in the x, y, z directions respectively 
x streamwise direction 


lateral direction 



z 


direction normal to the plate surface 


Greek Symbols 

S 0 hydrodynamic boundary layer thickness at the step 

y ratio of specific heats (isentropic exponent) 

p density 

p dynamic viscosity 

t w wall shear stress 

C frw2 ~ T wl )/(f w 2 — T aw ) 

Subscripts 

0 stagnation 

1 plate 

2 gauge 

A area-averaged 

aw adiabatic wall (recovery) 

L length-averaged 

x local 

w wall 

oo freestream 



1. Introduction 


Cylindrical flush mounted heat flux gauges (such as Schmidt-Boelter and plug type 
gauges) are frequently used on launch vehicles for measuring convective and radiative 
incident heat flux. The presence of the heat flux gauge necessarily alters the temperature and 
heat flux distribution on the gauge surface [1]. These gauges (normally metallic such 
aluminum, stainless steel, copper, etc.), when flush mounted in ablators or insulators of 
relatively low thermal conductivity (forming the vehicle surface), are subjected to wall 
temperatures considerably lower than that of the surrounding vehicle skin material (by 
several hundred degrees). On account of this surface temperature discontinuity (mismatch), 
heat exchange occurs between the gauge and the surrounding material by conduction, and 
the thermal boundary layer is altered (distorted), producing changes in the convective heat 
transfer coefficient from the fluid to the gauge surface. In convection, the thermal history of 
the boundary layer is carried with the fluid and affects the downstream heat transfer [1]. The 
changes in surface temperature and the departure of the heat transfer coefficient from the 
isothermal value demand that the heat flux measured by the calorimeter (gauge) deviates 
considerably from the true heat flux that would be measured if the gauge is formed of the 
same material as that of the insulation [2-3]. Experiments and analyses [4-10] have shown 
that a variable surface temperature distribution can produce a marked increase or decrease in 
the local and average convective heat transfer rates to a surface in laminar and turbulent flow. 

Correction factors have been proposed in the past to account for the departure of the 
measured heat flux from the true heat flux. Rubesin [4] and Reynolds et al. [5-6] presented 
two-dimensional integral solutions for the local heat transfer coefficient in incompressible 
flow past a flat plate with zero pressure gradient and a stepwise surface temperature with the 
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assumption of power law profiles for velocity and temperature. These solutions are based on 
the superposition of well-known unheated starting length solutions [2,11]. Recently, Mukerji 
et al. [12] extended the integral solution of Reynolds et al. [5,6] with an empirical correction 
(based on low speed data) to accommodate the effect of the viscous sublayer on the gauge 
surface (see Appendix A.3). Such a correction increases the influence of the temperature 
jump especially for small gauges, and offers a considerable improvement. On the basis of 
closed from integral solution due to Reynolds et al. [5,6], Westkaemper [13] derived a 
length-averaged heat transfer coefficient for the heat flux gauge. 

Many simplifying assumptions are considered in the development of the integrals. The 
l/7th power law solution is valid only far downstream of the step, as the profiles are 
developed at relatively large distances downstream of the step. Very near the temperature 
step, the thermal boundary layer is confined to the viscous sublayer. Artificial diffusivity in 
the sublayer (on account of l/7th power law profiles throughout), which is the most 
important part of the boundary layer, produces error in the region immediately downstream 
of the step in temperature [11]. Also, the boundary layer equations may not strictly apply 
near the leading edge of the discontinuity. Thus the accuracy of the integral solution at 
distances within a few boundary layer thicknesses is questionable. It is 
known that the integral solution due to Reynolds et al. [5,6] is known to underpredict 
experimental data by as much as 25 percent at jc/< 5 0 =1.0, and exceeds 30 percent for 

xl8 0 = 0.25 [12]. Another major deficiency of these integral solutions is connected with the 

assumption of constant thermal properties, which can lead to significant errors when the 
difference between the plate and the gauge surface temperatures becomes large, such as 


occurs at high Mach numbers encountered in flight conditions. In such circumstances 
solutions of the Navier-Stokes equations are necessary. 

Kandula and Reinarts [14] carried out a 2-D Navier-Stokes CFD analysis with variable 
properties to investigate thermal boundary layer corrections for the local heat transfer 
coefficient. It was shown that the correction factor for convective heat transfer coefficient 
increases with the wall temperature. To account for the circular geometry gauge, a quasi-two 
dimensional analysis was also carried out, wherein at any lateral plane of the gauge, the flow 
is assumed quasi-two dimensional, so that area-averaged heat transfer correction factor 
incorporating 3-D effects approximately can be estimated. Comparisons were then presented 
between the quais-2 D results from CFD (variable properties) and the constant property 
integral solutions, signifying the importance of property variations. 

The assumption of quasi-2D approximation can be in error as the size of the heat flux 
gauge is reduced in view of the preponderating effects of the three-dimensionality of the 
thermal boundary layer. In such circumstances, it is imperative to consider the full three- 
dimensional thermal boundary layer effect on the heat flux correction factors. To the authors’ 
knowledge, no such study has been reported previously. 

In this report, three-dimensional Navier-Stokes simulation of the thermal boundary 
layer has been carried out for the plate-gauge system subjected to a stepwise surface 
temperature discontinuity. The effect of various relevant parameters on the correction factors 
has been studied. Comparisons of the 3D analysis with quasi-2D CFD results and integral 
solutions have also been performed. The effect of three-dimensionality of the thermal 
boundary layer and of property variations on the convective heat flux corrections for the 
gauges has been investigated. 
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2. Three-Dimensional CFD Analysis 

2.1 Physical Assumptions 

The following physical assumptions are made in the CFD analysis. 

a. The heat flux gauge is in steady state (water-cooled). 

b. The gauge surface is idealized as isothermal. For example, a Schmidt-Boelter gauge is not 
strictly isothermal even when it is water cooled [15]. 

c. The external flow is uniform, and pressure gradient, surface curvature and separated flow 
are absent. 

d. The Prandtl number of the fluid is taken as 0.7. 

e. Both the plate and the gauge surfaces are smooth, with the gauge flush-mounted into the 
plate. 

2.2 Grid System 

The overlapping grid system for the plate/gauge configuration is shown in Fig. la, 
indicating the plate grid (118x72x84 size in the streamwise, lateral and normal directions) 
and the gauge grid (69x69x63 size in the radial, circumferential and normal directions). The 
total number of grid points is about 9.9xl0 5 . Individual grids for the plate and the gauge are 
generated. The inter-grid communication is provided by the Pegasus code, Benek et al. [16]. 
Fig. lb displays a side view of the grid. A plan view (partial) of the grid is provided in Fig. 
lc. 

An upstream inviscid plate length ( -0.5 < x/ L < 0 ) is considered. The plate leading edge 
and trailing edges are located at jc/L = 0and x/L = 2 respectively, with the gauge surface 
centered at x/L = l . The grid normal extent is taken as z / L = 0.9 . For flow resolution 
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purposes, the grid is clustered in the normal direction (near the wall) and in the axial 
direction (near the leading edge and near the surface temperature discontinuity). 

2.3 Flow Solution 

The flow solution (for density, velocity and temperature distribution) has been obtained 
by the OVERFLOW compressible Navier-Stokes CFD code [17,18]. A zonal two-equation 
k-co SST (Shear Stress Transport) turbulence model due to Menter [19] has been considered. 
The k-a> model has been validated for boundary layers and free shear layers [20]. 

Freestream boundary conditions are applied at the inflow boundary (x/L = - 0.5 ), top 
boundary and at the lateral boundary away from the symmetry plane. Lateral symmetry is 
used at the symmetry plane. An extrapolation boundary condition is considered at the 
outflow (xl L = 2). Viscous wall condition is prescribed at the plate and the gauge surface. 
Isothermal conditions are specified for the plate surface ( T = T W] ) and the gauge surface 
(T=T w2 ). 

The steady state solution is obtained by multi-gridding and local time-stepping. The 
algorithm used is diagonalized ADI (Alternating Direction Implicit), with central 
differencing for fluxes. Sutherland correlations are considered for the temperature 
dependence of viscosity and thermal conductivity of air [21, 22]. 

j 

2.4 Flow and Geometry Parameters 

Solutions were obtained for a freestream Reynolds number of 1E6, and freestream 
Mach numbers of 2 and 4. Values of the gauge radius to plate length ratio R/L= 0.005, 0.01, 
0.025, 0.05, 0.075 and 0.1 are considered. The freestream temperature is 288 K (519 R), and 
the gauge surface temperature is held at 333 K (600 R). The plate wall temperature is varied 


8 



in the range 389 K to 1944 K (700 R to 3500 R), and comprises values on either side of the 
adiabatic wall (recovery) temperature. The recovery temperature is computed from 

Z™L = i + r IziM 2 (1) 

T„ 2 

where the turbulent recovery factor is obtained from 

r = Pr m (2) 

With an assumed value of Pr = 0.7 , the above formula yields a recovery temperature of 494 K 
for M = 2, and llllKforA/=4. 

3. Results and Comparison 

As indicated in the assumptions, all the results presented here with regard to CFD and 
the integral models are limited to uniform flow past a flat plate, where pressure gradient, 
surface curvature or separated flow are absent. The computations are performed on an 8- 
processor (R- 10000) SGI Origin-2000 machine. 

3.1. Solution Convergence 

Fig. 2 shows the convergence history for the solution residuals in a typical case. This 
residual history serves as a qualitative measure of the convergence to steady state, and 
suggests that convergence in residuals (several orders of magnitude drop in residuals) for the 
two grids is approached in about 2000 time iterations. 

3.2 Local Stanton number Distribution 

The local Stanton number distribution in the symmetry plane for a particular case is 
illustrated in Fig. 3a for Rl L = 0.01 mdM = 4 , T wl IT m = 0.5 , and T w2 !T m = 0.3 . A jump in 

the local Stanton number on the gauge surface (0.99 <jc/l< 1.01) is noted, as is to be 
expected. Downstream of the gauge, the recovery in the local Stanton number is established. 
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The close agreement between the plate solution and the gauge solution in the overlapping 
region suggests that the overlapping grid scheme is satisfactory. 

A schematic of the thermal boundary layer in the vicinity of the discontinuity for the 
symmetry plane is illustrated in Fig. 3b. The emergence of a new thermal boundary layer past 
the discontinuity is illustrated. Fig. 3c highlights the Stanton number distribution on the 
gauge surface. The extent of the three-dimensionality of the thermal boundary layer is 
evident near the leading edge and the trailing edge of the gauge. 

3.3 Heat flux correction factors 

The effect of variation of properties on the heat flux correction factors is illustrated in 
Fig. 4a for the case of M = 4 and R/ L = 0.01 . The constant property quasi-2D solutions for 
the integral method and the modified integral method [12] are compared with the quasi-2D 
CFD results for both constant and variable properties, and with the 3D CFD variable property 
solution. In the legend, CP stands for constant properties, and VP for variable properties. In 
general, a linear variation of q 2 / q x with 0 is noted, as predicated by the integral result of Eq. 
[A.5]. It is seen that the modified integral solution shows better agreement with the 2D CFD 
constant property solution. The variable property 2D CFD solution falls between the constant 
property 2D CFD solution and the modified integral solution. The 2D CFD solutions reveal 
that as the plate temperature to recovery temperature difference increases, the effect of the 
property variation becomes more pronounced. 

Fig. 4a also suggests that the 3D CFD solution for the heat flux correction factor 
exceeds that obtained by quasi-2D CFD solution. The effect of three dimensionality is seen 
to increase with an absolute increase in 0 . The physical mechanism for the enhailced heat 
flux in three dimensions may be explained as follows. As the chord length of the gauge 


decreases away from the center-line, the quasi-2D assumption becomes more and more 
unrealistic, and the leading edge effect increases the overall 3D heat flux in the presence of a 
step. Data obtained by Mukeiji et al. [12] under low turbulence (4 percent) and high 
turbulence (7-10 percent) freestream turbulence level suggest that spanwise transport of heat 
through turbulent diffusion is not an important effect for a 2-D heated element (increases the 
heat transfer rate by at most 5 percent). 

Fig. 4b shows the corrections factor variation as a function of T wl /T aw for T wl IT aw < l . 

The results reveal that the deviation between 2D and 3D CFD increases as the wall 
temperature increases, indicating the importance of three-dimensional thermal boundary 
layer. The modified integral solution is seen to offer a considerable improvement over the 
integral solution when compared with the CFD results. 

Comparisons for the correction factors for M = 0.2 and R / L = 0.01 are presented in Fig. 

5. The general trends are seen to be similar to those expressed in Fig. 4a. A direct 
comparison between 3D CFD results at M =4 and M = 2 is displayed in Fig. 6, showing the 
effect of Mach number on the heat flux correction factors. The results reveal that the 
correction factor increases as the Mach number increases, as is to be expected owing to the 
role of compressibility and property variations. It is expected that the agreement would be 
closer if the values of T w2 IT aw were kept constant in both cases. 

Calculations for M = 4 at R/ L = 0.1 (Fig. 7) indicate that the heat flux correction factors 
are smaller than those for R/ L = 0.01 . Although the trends are similar those for R / L = 0.01 , it is 
observed that there is no appreciable difference between the 3D and the 2D CFD results, 
suggesting that the effect of three dimensionality becomes negligible for R/ L = 0.1 . 
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The heat flux correction factors for M = 4 at R / L = 0.005 are depicted in Fig. 8. 

AtR/ L = 0.005 , the heat flux corrections factors are found to be higher than those 
for R/L = 0 . 01 , as is to be expected. Again the modified integral solution agrees better with the 
2D CFD solution as compared with the original integral solution. It is remarkable that the 3D 
CFD results are smaller than those from 2D values in the range of T wl n aw from 0.3 to 0.45. 
The physical basis for this behavior is not clear. 

3.4 Comparison with test data 

In an effort to validate the predictions, comparisons are made with limited test data 
recently reported by Strobel et al. [23] at M = 4 and Re= 1.2E6. The data were obtained from 
ground tests at NASA Marshall hot gas facility (HGF) employing hydrogen-air combustion. 
Heat flux data were obtained from Schmidt-Boelter gauges (4.76 mm diameter, aluminum) 
supplied by AEDC (Arnold Engineering Development Center), and reference thin-skin 
calorimeters (2.54 mm thick and 1 14 mm diameter; 17-4 PH stainless steel) with the latter 
designed to eliminate 2D conduction. A range of heating conditions can be produced by 
varying combustor pressure and temperature. The test section is comprised of a flat plate 
(12.7 mm thick, 304.8 mm wide, and 501.7 mm long). Both bare steel plate and steel plate 
covered with acusil TPS (Thermal Protection System) insulating material are considered. The 

t 

heat flux gauges are not water-cooled. Instead, the plate-gauge system is operated in 
dynamic (transient) mode. The measured surface temperatures suggest that the temperature- 
time history becomes quasi-steady in a few seconds (the test duration is about 20 sec). 

The thin skin heat flux is found to be comparable regardless of the surrounding material 
(steel or acusil). There are small discrepancies between heat flux gauge output and thin skin 
response when gauges are mounted in stainless steel. However large discrepancies are noted 
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when the gauges are mounted in acusil on account of changes in the thermal boundary layer. 
The original data are reported in terms of cold wall heat flux (T w = 300 K), and are converted 
to hot wall heat flux. Considering uncertainties in surface temperature measurement, material 
properties, and dimensions, and errors associated with lateral heat conduction loss, the 
measured heat flux data are believed to be within about seven percent accuracy. 

Fig. 9 shows the variation of q 2 1 q x as a function of R/ Las predicted by various 
theories and their comparison with test data. The data at 10.2 atm stagnation pressure 
correspond to wall temperatures of T wl /T aw =0.508 and T^/T^ =0.303, and those at 13.93 atm 
stagnation pressure correspond to T wl / T aw = 0.515 and T w2 lT aw = 0.296 . The predictions suggest 
that q 2 / q x decreases with an increase in the value of R/L. The 3D CFD results are higher 
than 2D CFD values, and the deviation increases with a decrease in R/L. The test data 
corresponding to R/ L= 0.01 and 0.025 are seen to lie between the 2D CFD and 3D CFD 
results. In general the CFD results exceed considerably the constant property integral 
solution. For the range of conditions considered here, the integral solution of [5,6] is within 
5-25 percent lower than the CFD results. The modified integral solution [12] shows better 
agreement with the data and the CFD solutions, being within about 10 percent lower than the 
CFD. These comparisons thus tend to highlight the importance of the three dimensional 
boundary layer effects. 

Referring to the accuracy of the turbulence model employed here, the two-equation 
k - (O model has been widely incorporated in 3D CFD codes and validated for complex three- 
dimensional configurations (for example, see Refs. 28 and 29). Calculations by the authors for 
local Stanton number for an isothermal flat plate at Re=10 6 and M = 0.2, the model is accurate 
within 2 to 6 percent of standard correlations over a wide range of temperatures. Thus for 



boundary layer type flows considered here (where pressure gradients and flow separation are 
absent), it is estimated that the uncertainty in the heat flux correction factor due to turbulence 
model is about 6 percent. 

4. Conclusion 

The three dimensional Navier-Stokes solution of the heat flux gauge subjected to a 
surface temperature discontinuity has led to some insight into the uncertainty in heat flux 
measurement. It has been shown that the heat flux correction increases with an increase in 
plate temperature relative to the gauge temperature. It also increases with decreasing value of 
the gauge radius, and is relatively weakly dependent on the flow Mach number. In general, 
the effect of variable properties becomes significant with an increase in the temperature 
differential between the plate and the gauge. The results suggest that the effect of three- 
dimensionality tends to increase the heat flux correction considerably above the integral 
predictions. The modified integral correlation accounting for the sublayer effect performs 
much better than the integral solution. For relatively large values of R/L (say R/L = O.lor 
above), the three dimensional thermal boundary layer effects become insignificant. For 
relatively small values of R/ L (=0.005), there is seen to be a crossover point between the 3D 
CFD and 2D CFD results. 

Comparisons with limited test data available agree favorably with the three dimensional 
thermal boundary layer results, and confirm the importance of three dimensionality of the 
heat flux gauge and of property variations. Additional comparisons are recommended to 
further validate the accuracy of the model over a wider range of test conditions. As the heat 
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flux corrections are generally huge, it can be concluded that the water-cooled heat flux 
gauges should never be considered for measurement of convective heat flux. 
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Appendix A. Two-dimensional integral solutions 

A.l Local heat transfer coefficient 

Consider the flow past a two-dimensional flat plate with a surface temperature 
discontinuity 

T w =T W i 0 <x<L 

= T W 2 x>L 


(A.l) 


With the aid of an integral method with assumed power law profiles for velocity and 
temperature, Rubesin [4] and later Reynolds, Kays and Klein [5,6] obtained an expression for 
the local heat transfer coefficient h(x,L) in turbulent incompressible flow in the form 


^-4 + < 

h(x, 0) 


fll- 






ffl2 


x > L 


where b = {t w1 -T aw )/(T w2 -T aw \ ( = {t w2 -T wl )/( t w2 -T aw ) 

Here, h(x,o) denotes the local heat transfer coefficient on an isothermal plate with constant 
wall temperature T wl . The exponents m, and m 2 are 39/40 and -7/39 respectively from Rubesin 
[4], and 9/10 and -1/9 from Reynolds [5, 6], which is valid over a wider range of Reynolds 
number. For a single step in temperature (T wl = T„,T w2 = T w ), the solution due to Reynolds et 


(A.2) 

(A.3) 


al. [5,6] is expressed by 


St x = 0.0287 Re; 0 2 Pr ^ 4 [l - (£7 x) 9/10 ]~' /9 


(A. 4) 


The corresponding local heat flux ratio is expressed by [9] 


(A.5) 


where 

and 


q i = M*. i-Xr.1 -T„l q,= h(x,0jT wl -T„) 
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For very small values of L/x, q 2 l , which depends on temperatures only, and 
independent of geometric parameters [9]. 


A. 2 Length averaged heat transfer coefficient 

Based on the integral solution of [5,6], Westkaemper [12] derived a length-averaged 
heat transfer coefficient h L over the heat flux gauge as: 


V L =- 


hL 


W 




= bF(L/w)+£H{L/w) 


where h L =h(w,L) is given by 


(A.6) 


h L = 


' w w 

. Q=\ q{x)dx = \ h{x, l){T w2 -T aw )dx 




The factors F and H are defined by 


5 [i-(l/w) 4/5 


F = c k 


4 (l-(L/W)) 


H=Ck f flT/w) ^ /L)9/ '° -J ' 9 -Hl/w) 

where c k = l . Later, Knox [24] pointed out an error in Westkaemper’s equations for F 


and H . Considering that 


fW + L ^ 


,0 


V ^ J 




A + L/W j 


(A.7) 


(A. 8) 


it was shown that the correction is provided by 

c* =[2/(l + L/W')] 1/5 


(A.9) 
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A. 3 Modified Integral Solution due to Mukerji et al. [12] 

Based on STA7 2D boundary layer finite difference code in conjunction with low speed 
data, Mukeiji et al. [12] recently proposed the following correlation, as an extension of the 
integral solution of Reynolds et a. [5,6]: 


5r(x) = 0.0287 Re^ 2 Pr^ 4 - 

1- 

^y /i0 ' 

-1/9 

+ (p(^,x,Rt^)> 


L 

<x) 



(A. 10) 


where 


</> = (- 0.0139 ln(Re^ )+ 0.246)f ^ 


^/l -(£/*) 


and represents a correction to the Reynolds correlation. For large values of x/S 0 (where 
^represents the hydrodynamic boundary layer thickness at the location of the temperature 
step £, and xthe downstream distance from the step), this relation reduces to Eq. [A.4] due 
to Reynolds et al. [5,6]. The square root term controls the shape of the curve near the step. 
The correlation for local Stanton number agrees within 5 percent of STAN7 for x /S 0 = 0.05 

to 10, and Re^=5xl0 5 to 3xl0 6 . It predicts low turbulence intensity (0.4 percent) rectangular 
calorimetric data within two percent. 

It may be pertinent to point out that the near-step region has been investigated in [25,26 ] 
for a step change in surface temperature, and in [27] for a step change in surface heat flux. 


Appendix B. Quasi-two-dimensional extension 

Kandula and Reinarts [13] considered a quasi-two dimensional extension to 
approximately account for the cylindrical geometry of the heat flux gauge. At any lateral 
plane of the gauge, the flow is assumed quasi-two dimensional, so that an area-averaged heat 
flux correction factor h A incorporating the three-dimensional effects can be estimated based 
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on the two-dimensional results (both variable property CFD and constant property integral 
solutions): 


where 


V A 


-v*f- 





R 

fr L {W',L'X2x')dy', 
y '= o 


x = R cos0, y = Rs\nd 


(B.l) 

(B.2) 


Calculations suggest that the quantities h A and h L denoting quasi-2D and 2-D corrections 
do not appreciably differ from one another. For instance, for R/ L = 0.01 at 
M =4, R = 1E6, T wl /T aw = 1.6, the quantity h A is only 2.5% higher than h L . 
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Flow 



Plate leading edge 


Fig. la. Overall view of the 3-dimensional grid system. 
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Plate wall at T w i 
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Plate leading edge 


Gauge at T W 2 


Fig. lb. Side view of the grid system. 



Solution residual 



Fig. lc. Plan view of the grid system in the neighborhood of the gauge. 



Fig. 2 Convergence history of solution residuals. 
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Fig. 3a. Local number Stanton number distribution in the symmetry plane. 



Fig. 3b Thermal boundary layer alteration near the temperature discontinuity 
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Fig. 3c Stanton number distribution on the gauge surface. 



<t> 

Fig. 4a Comparison of heat flux correction factor for M=4 and R/L=0.01. 
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Fig. 4b Comparison of heat flux correction factor for M=4 and R/L=0.01 (contd.). 
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Fig. 5 Comparison of heat flux correction factor for M=2 and R/L=0.01. 
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Fig. 6 Comparison of heat flux correction factor from 3D CFD for M=2 and 4 
At R/L=0.01. 
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Fig. 7 Comparison of heat flux correction factor for M= 4 at R/L=0.1. 
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Fig. 8 Comparison of heat flux correction factor for M= 4 at R/L=0.005. 



Fig. 9 Comparison of heat flux correction factors with test data for M= 4. 
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